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Spectral Decline of PeV Neutrinos from Starburst Galaxies 
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Starburst galaxies represent one of the most plausible origins of the cosmic high-energy neutrino 
flux recently discovered by IceCube. At ~ PeV energies, the neutrino flux from starburst galaxies 
is expected to exhibit a characteristic spectral break due to cosmic-rays escaping the galaxy. We 
examine the ’smearing’ of this spectral break by a population of starburst galaxies with varying 
properties. We incorporate galaxy distribution w.r.t. star-formation rate and redshift. Our results 
(i) show characteristic spectral softening in IceCube’s energy band; (ii) resolve the conflicting obser¬ 
vations of soft neutrino spectrum and diffuse gamma-ray flux observed by Fermi-LAT; (iii) constrain 
the properties of the magnetic fields in starburst galaxies. 


I. INTRODUCTION 

IceCube has recently discovered a quasi-diffuse flux of 
astrophysical neutrinos [1-3]. The origin of these neutri¬ 
nos is currently unknown. Unraveling the source will help 
understand the origin of cosmic rays, the mechanisms be¬ 
hind particle acceleration, and the broader environment 
of the acceleration site. 

The astrophysical neutrino flux has been observed in 
the 20TeV - 3 PeV range to be ~ 10“® GeVcm“^s“^sr“^ 
per neutrino flavor [2, 4]. The spectrum of the flux is, 
however, currently uncertain. Observations seem to be 
inconsistent with dN,y/dEi, oc E~^, with T = 2 spectrum 
typically expected from Fermi acceleration [2, 4]. Con¬ 
straints using muon neutrinos from the Northern hemi¬ 
sphere detected by IceCube between 2010-2012 in the en¬ 
ergy range 330TeV-1.4PeV indicate an effective spectral 
index of Pefi « 2.2 [3]. Using a combination of IceCube 
searches with neutrinos observed in the 25TeV-2.8PeV 
energy band, Aartsen et al. [2] finds a best fit spectral 
index of Peff = 2.50 ± 0.09. Results are also compatible 
with a harder spectrum with a spectral break. Such spec¬ 
tral features are plausible for a variety of source types, 
including starburst galaxies [5, 6] or merger shocks in 
clusters of galaxies [7, 8]. 

Starburst galaxies represent one of the prime candi¬ 
dates for the origin of high-energy neutrinos [5, 9-13]. 
High star formation rate (SFR) is connected to the occur¬ 
rence of phenomena such as gamma-ray bursts [14] and 
supernova remnants [15] which are plausibly connected 
to the acceleration of cosmic rays [16, 17]. Even more 
importantly, strong magnetic fields and high interstel¬ 
lar gas densities in starburst galaxies make them largely 
opaque to most cosmic rays [10-12]. Cosmic rays pro¬ 
duced within starburst galaxies lose their energy through 
nuclear collisions, producing high-energy neutrinos. Only 
the highest energy cosmic rays can escape since they are 
less affected by magnetic fields in the galaxies. 
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FIG. 1. Expected quasi-diffuse neutrino spectrum from star- 
burst galaxies, adopting injection index F = 2. Results are 
shown for three different magnetic held models (see legend), 
Rmin being our benchmark model. Also shown is an empir¬ 
ical ht onto the benchmark spectrum (dotted line), and a 
power law ht onto the lOTeV-lOPeV energy interval (black), 
which is the energy band of astrophysical neutrinos detected 
by IceCube. Additionally shown are the observed IceCube 
and Fermi quasi-diffuse huxes, as well as the 7 hux expected 
given the benchmark neutrino spectrum. The inserted plot 
separately shows the benchmark spectrum and power-law ht. 


An additional constraint on the neutrino spectrum 
comes from the quasi-diffuse GeV gamma-ray flux ob¬ 
served by Fermi-LAT [18]. If cosmic neutrinos are pro¬ 
duced in hadronuclear interactions, as expected for the 
case of starburst galaxies and many other models, there 
is a corresponding gamma-ray emission that represents 
an independent probe of the neutrino spectrum. With 
the latest IceCube observations, the observed Fermi- 
LAT quasi-diffuse gamma-ray flux constrains the neu¬ 
trino spectral index toF < 2.1 — 2.2, given that the 
neutrinos are produced in pp interactions [19]. This con¬ 
straint holds even for a broken power law spectrum. To¬ 
gether with the current best fits on IceCube observations, 
requires either (i) at least partly photo-hadronic neutrino 
production or (ii) deviation from single power law spec- 
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trum. 

In this paper we study the neutrino spectrum ex¬ 
pected from starburst galaxies, accounting for their vary¬ 
ing properties and cosmic distribution. Our goal is to 
gain a more detailed picture of the spectral decline due 
the cosmic ray escape energy, which varies between star- 
burst galaxies. After reviewing the expected neutrino 
spectrum from a single starburst galaxy, we use SFR to 
recover corresponding typical values for galaxy proper¬ 
ties. We then consider the cosmic distribution of star- 
burst galaxies to derive the cumulative neutrino spec¬ 
trum. We finally compare the recovered deviations from 
the power law spectrum expected from the cosmic ray 
sources within galaxies with neutrino and gamma-ray ob¬ 
servations. 


II. NEUTRINO SPECTRAL CUTOFF IN A 
SINGLE STARBURST GALAXY 


We define the neutrino spectrum expected from a sin¬ 
gle starburst galaxy following [10]. Below a threshold 
energy Ath, the spectrum of cosmic rays in the star- 
burst galaxy is determined by their production spectrum, 
which we assume to follow a power law with spectral in¬ 
dex p, and the confinement time Tconf ■ In the Milky Way, 


(iogIv) <'> 

for cosmic ray energy Ap, with s « 0.6 [10]. 

The threshold i?th is determined by the relation be¬ 
tween the confinement time and the energy loss time 
of protons, assuming that the starburst lifetime is suf¬ 
ficiently long [10]. The energy loss time for hadronuclear 
interactions is rioss ~ (O.SruTppc)”^, where n is the inter¬ 
stellar nucleon density, and the inelastic nuclear collision 
cross section CTpp « 50 mb for the energy range of interest. 
The factor 0.5 is the collision’s inelasticity. Expressing 
Toss with gas surface mass density Eg m-puh, h being 
the galactic disk height, we have 


”7"loss 


= 7xio^( 
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where ha ~ 300 pc is the Galactic disk height. For the 
confinement time, we start with Eq. 1, and assume that 
it depends on the cosmic ray energy only through the 
Larmor radius (oc E/B). We further assume that the 
confinement time scales as the square of the galactic disk 
height since this is a diffusion process. We arrive at 


'^conf 


10^ 


Er, 


lOGeV 




yr 


( 3 ) 


where Bq 6 pG is the Galactic magnetic field strength. 

We can further utilize the connection between B and 
Eg. Primarily, B is estimated using the observed syn¬ 
chrotron emission of cosmic-ray electrons, although there 


are indications that the magnetic fields in starburst 
galaxies can be significantly stronger than this estimate 
[21, 22]. We follow McBride et al. [22] and denote this 
magnetic field estimate based on synchrotron emission 
as Bmin, since this represents a minimum strength. This 
will be our benchmark model. Based on McBride et al. 
[22], Bmin can be expressed as 


Bn,in = 10 
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Magnetic field strengths can also be derived from the 
radio luminosities of supernova remnants. We will de¬ 
note this estimate with Bgnr- Another possibility is that 
the magnetic energy density is equal to the pressure of 
a self gravitating disk in hydrostatic equilibrium, as a 
function of gas surface density. We will denote the esti¬ 
mate based on this connection with Bhyd- For these two 
other possibilities the conversions according to McBride 
et al. [22] are Bgnr = 6 x 10“"‘(Eg/gcm“^)°-®® G and 
Bhyd = 62 X 10-3(Eg/gcm-2) G. 

Using these conversions, we can now determine Bth 
by finding the energy for which tioss = Tonf- For our 
benchmark model with Bmin, we find 

=6.5x 10“ (^)” ’ TeV, (5) 


For the other two models, we find Eth(Bsnr) = 4 x 
103(Eg/gcm“2)2-2(/i//i(^)i-7TeV and Bth(Bdy„) = 6 x 
10^(Eg/gcm“^)^-^(/i//iG)^'^TeV. Note that Eth is the 
threshold for cosmic rays; the corresponding threshold 
in the neutrino spectrum is ~ O.ObBth since neutrinos 
are produced with about 5% of the proton energy in in¬ 
teractions. 


III. STARBURST GALAXY CHARACTERISTIC 
PARAMETERS FROM STAR-FORMATION 
RATE 

We saw above that the spectral features of a starburst 
galaxy can essentially be expressed with one parameter: 
Eg. The weight of each galaxy in an ensemble spectrum, 
however, also needs to be taken into account. Further, 
there is no readily available catalog or model for the cos¬ 
mic distribution of Eg for starburst galaxies. To over¬ 
come this, here we derive a connection between SFR and 
Eg. SFR can be used as a measure of the cosmic-ray emis¬ 
sion from a starburst galaxy. It is also more straightfor¬ 
ward to observationally determine than Eg. The cosmic 
distribution of starburst galaxies w.r.t. SFR has been 
characterized [23]. 

Starting with the SFR of a starburst galaxy (deter¬ 
mined, e.g., from its far infrared emission [24]), we first 
convert SFR to SFR surface density Esfr- For galac¬ 
tic half-light (or effective) radius B, we define Esfr = 
SFR/7rR2 (e.g., [25]). 
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To determine the characteristic value of R, we turn to 
the Schmidt relation [26] 

SsFR = 30 -^ ) Mq yr"^ kpc"^ (6) 

Vgcm-^y 

Substituting Eg = M^j-KR? into Eq. 6 with Mg being 
the galaxy’s molecular gas mass, we can express R as 


= 2 X 10"^® 


SsFR 


Mq yr“i kpc 


- 1.2 




The characteristic radius can also be used to determine 
the characteristic disk height. Following Law et al. [27], 
we adopt h = 0.3R. 

We can further derive the typical value of Mg as 
a function of SFR. Assuming typical dynamical time 
T'dyn ~ 10^ yr, the conversion based on Genzel et al. [25] 
(their Eq. 8) is 


Mg = 5 X 10 ^ 


SFR 

Mgyr-i 


0.73 

Mq. 


( 8 ) 


Eqs. 6, 7 and 8 together express the typical value of 
Eg corresponding to a given SFR for starburst galax¬ 
ies, requiring no other parameter. Note that the above 
equations represent average relations between parame¬ 
ters. For an example, we take starburst galaxy Arp 220, 
whose gas density Eg^P^^^ = 5.3gcm“^ is quoted in [23]. 
The SFR of Arp 220 is ~ lOOMQyr”^ [28], for which our 
conversion gives Eg « 14gcm“^, a good estimate given 
the uncertainties of the measured properties and the vari¬ 
ations between starburst galaxies. 


IV. STARBURST GALAXY POPULATION 


With the above conversion from SFR to Eg and from 
Eg to Eth, we are now able to determine the expected 
neutrino spectrum for a starburst galaxy with known 
SFR. We now look at the distribution of starburst galaxy 
SFR values that allows the derivation of the cumulative 
neutrino spectrum. 

We adopt the starburst galaxy infrared luminosity 
function by Sargent et al. [23]. Sargent et al. find 
that the luminosity distribution is well described by the 
Schechter function. We convert their infrared luminos¬ 
ity values at redshift z ^ 1.1 (see Table 1 of [23]) to SFR 
based in Kennicutt [24] , and fit the derived SFR distribu¬ 
tion 4 >sfr with a Schechter function with fixed index a = 
1.4 (see [23]), finding parameters 4>* = 2.2 x 10® Mpc“^ 
and L* = 625MQyr“^. 

The SFR distribution also changes with redshift. Fol¬ 
lowing Sargent et al. [23], we take to be constant out 
to z = 1, and for z > 1 we adopt scaling with (l-|-z)“^-^°. 
We scale L* with (1 -|- z)^ ®. 

V. RESULTS 


A. Overall neutrino spectrum at Earth 


The SFR distribution obtained above is converted to 
the distribution 4>Eth of cosmic ray Fth weighted with 
galaxy SFR. This latter incorporates the assumption that 
cosmic ray production is proportional to SFR. Integrat¬ 
ing over SFR and redshift gives 


^Eth(A(h) oc 



/ ‘J’SFR 
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where $sfr = $sfr(SFR, z) and Fth = Ath(SFR). $Eth 
is then converted to neutrino spectrum 


dN,(E,) 

dEi, 


( 10 ) 


The resulting spectrum is shown in Fig. 1, for all three 
magnetic field models considered above. We see that the 
spectrum deviates from the single power law distribution 
of cosmic ray acceleration sites {E~^) within the sensitive 
energy band of IceCube. To describe this softening, we 
find the empirical fit 


dE,^ 


oc Aj, ^ exp 


E, 

UPeV 


0 . 43 ' 


( 11 ) 


for our benchmark model (see Fig. 1). The deviation 
appears to be slower than exponential. Our results are 


consistent with that of [5] who also considered the neu¬ 
trino spectrum from a population of galaxies, albeit with 
somewhat different underlying assumptions. 


B. Comparison to IceCube spectral fits 

We see that the spectral softening occurs in the Ice- 
Cube’s sensitive energy band for the benchmark mag¬ 
netic field model. This is in agreement with the predic¬ 
tion of Loeb and Waxman [10], who came to this conclu¬ 
sion assuming a population of identical starburst galaxies 
with characteristic properties (using the same benchmark 
model for magnetic fields). For the alternative magnetic 
field models, the spectral softening is above IceCube’s 
band. 

If a single power law is fit on the observed distribution 
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in IceCube’s sensitive band, the observed Fobs will be 
greater than expected for a single galaxy below the cutoff 
energy. Fitting a power law on the neutrino energy band 
lOTeV-lOPeV for the benchmark model, we find Ffit ~ 
2.2. This spectral fit depends on the energy boundaries 
considered for the fit. 


C. Comparison to GeV gamma-ray diffuse flux by 
Fermi-LAT 

If neutrinos are generated in pp interactions, there will 
be a corresponding gamma-ray emission as well. We cal¬ 
culated the expected gamma-ray flux corresponding to 
our benchmark neutrino flux similarly to the calculations 
of Murase et al. [19], taking into account the gamma-ray 
spectral attenuation due to propagation over cosmic dis¬ 
tances. Results are shown in Fig. 1, in comparison with 
the diffuse gamma-ray flux detected by Fermi-LAT [29]. 
We find that the obtained gamma-ray flux is compati¬ 
ble with the observed diffuse gamma-ray flux, with the 
gamma-ray flux being ^ 60-100% of the observed flux 
around ^ 100 GeV. 


D. Constraints on magnetic fields in starburst 
galaxies 

We see in Fig. 1 that spectral softening occurs within 
the energy band 10 TeV-10 PeV for the benchmark model 
B = Rmin- For the other two models, the neutrino spec¬ 
trum essentially follows a power law distribution within 
lOTeV-lOPeV, and softens only at higher energies. Us¬ 
ing the empirical fitting function in Eq. 11, we find that 
the characteristic energy of spectral softening is 14 PeV 
for B = Rmin, 110 PeV for B = i?snr, and 100 PeV 
for B = i?hyd- The reconstructed softer spectrum by 
IceCube therefore supports B = Rmin- 


VI. CONCLUSION 

We characterized the expected high-energy neutrino 
spectrum from all starburst galaxies by incorporating the 
galaxies’ luminosity function and its cosmic evolution. 
To accomplish this, we used starburst galaxies’ infrared 
luminosity to derive characteristic properties that deter¬ 
mine the neutrino spectrum and luminosity of individual 
galaxies, such as star formation rate and the magnetic 
field strength of the interstellar medium. Are conclusions 
are the following. 

• Spectral decline: Due to the distribution of star- 
burst properties, the expected neutrino spectrum soft¬ 
ens within the IceCube energy band 10 TeV-10 PeV) 
compared to a single power law spectrum, but the de¬ 
cline is not as sharp as the exponential cutoff expected 
for a single starburst galaxy. We find the empirical fit 


dN^/dE oc E ^ exp[(i?/£'o)^] with f3 « 0.43, represent¬ 
ing a slower-than-exponential break. 

• IceCube spectral fit: This has important implications 
to the spectral fit of IceCube observations. Namely, the 
observed spectral index within ~ 10 TeV-10 PeV range 
will be softer, Fes ^ 2.2, compared to the injection in¬ 
dex F = 2 for a single starburst galaxy below the cutoff 
energy. This observed softening is consistent with the 
spectral slope obtained by IceCube above ^ 100 TeV, al¬ 
though there may be additional components in the 10-100 
TeV range [20]. 

• Fermi-LAT constraints: The expected diffuse gamma- 
ray background at Earth corresponding to the neutrinos 
flux is compatible with the spectral constraints from the 
GeV gamma-ray background observed by Fermi-LAT. 
The derived gamma-ray flux is ~ 60-100% of the observed 
flux. Although there is significant tension with the stan¬ 
dard blazar interpretation of the diffuse gamma-ray back¬ 
ground, the fact that the two diffuse fluxes are compatible 
may further corroborate starburst galaxies as the poten¬ 
tial origin of the observed cosmic high-energy neutrino 
flux. Since cosmic rays are likely injected into starburst 
galaxies by energetic transient sources such as gamma- 
ray bursts and supernovae, this scenario can be promising 
for multimessenger transient observations, such as those 
using gravitational waves or gamma rays [30, 31]. 

• Magnetic field constraints: We investigated alternative 
magnetic field models for starburst galaxies with respect 
to their effect on the expected neutrino spectrum. We 
find that for the benchmark magnetic field based on the 
observed synchrotron emission from cosmic ray photons 
predicts neutrino spectral softening in the energy band of 
recent IceCube results, while other magnetic field models 
predict significantly higher neutrino break energies, mak¬ 
ing the spectrum in the IceCube band essentially an un¬ 
broken power law. The confirmation of spectral softening 
in the ^ 10 TeV-10 PeV band compared to at TeV ener¬ 
gies by IceCube can therefore provide useful constraints 
on the interstellar magnetic field strength of starburst 
galaxies. 

We foresee two essential future improvements of our 
calculations. First, there can be more direct measure¬ 
ments of starburst galaxy properties, such as magnetic 
fields, for a large number of starburst galaxies, which can 
provide a more accurate estimate than the derived quan¬ 
tities that we used above. Second, our simple starburst- 
galaxy model can be improved beyond having a sharp 
spectral cutoff due to the escape of cosmic rays by a more 
realistic accounting for the fraction of cosmic rays escap¬ 
ing as a function of energy. These improvements can 
then be directly incorporated in the rest of the analysis 
presented above. Finally, the general concept presented 
here, namely the spread of a spectral feature due to a 
population of diverse sources, can be similarly applied to 
source candidates other than starburst galaxies. 
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